APS/123-QED 



Enhanced quasiparticle heat conduction of the multigap superconductor Lu 2 Fe 3 Si5 
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The thermal transport measurements have been made on the Fe-based superconductor Lu2Fe3Sis 
(T c ~ 6 K) down to a very low temperature T c /120. The field and temperature dependences of the 
thermal conductivity confirm the multigap superconductivity with fully opened gaps on the whole 
Fermi surfaces. In comparison to MgB2 as a typical example of the multigap superconductor in a 
p-electron system, Lu2Fe3Sis reveals a remarkably enhanced quasiparticle heat conduction in the 
mixed state. The results can be interpreted as a consequence of the electronic correlations derived 
from Fe 3d-electrons. 
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Multigap superconductivity (MGSC) is the existence 
of superconducting gaps with significantly different mag- 
nitude on distinct Fermi surfaces. This phenomena 
has been realized in wide variety of materials including 
MgB 2 [J, NbSe 2 and the heavy fermions 0-d, point- 
ing to its universality underlying the superconductivity. 
One consequence of the MGSC is the ability to excite 
low-energy quasiparticles (QPs) due to the presence of 
the small gap, providing unusual features in the mixed 
state. Even though the MGSC has been extensively stud- 
ied so far, the purpose of the study has been mainly 
focused on the explanation for the anomalous proper- 
ties observed in the superconducting state. Recent dis- 
covery of the iron-pnictide superconductors Q, however, 
have offered further insight into the MGSC especially in 
the correlated electron system because the electronic in- 
teractions and the multiband structure are essential for 
the pairing mechanism in this system Q. In that sense, 
a question of how the electronic correlations affect on 
the MGSC is a fascinating issue to be addressed, which 
could not be examined on MgB 2 . Unfortunately, lack of 
high-quality stoichiometric samples and/or the high up- 
per critical held prevent detailed studies of the MGSC in 
the iron-pnictides 0j- 

Lu2Fe3Si5, a ternary-iron-silicide, is an another ex- 
ample of the Fe-based multigap superconductor with 
T c ~ 6 K [8], which crystallizes in the tetragonal struc- 
ture consisting of a quasi-one-dimensional iron chain 
along the c axis and quasi-two-dimensional iron squares 
parallel to the basal plane Band structure calcula- 
tions predict that Fermi surfaces consist of two holelike 
bands and one electronlike band, and each band has a 
contribution from the Fe 3<i-electrons In the two 
holelike bands, some parts of the Fermi surfaces have 
the different dimensionality reflecting the crystal struc- 
ture Q. Importantly, the Fe 3<i-electrons are responsi- 
ble for the superconductivity of Lu 2 Fe 3 Si5, as suggested 
by the absence of superconductivity in the isoelectronic 
Lu 2 Ru3Si5 and Lu 2 0s3Sis jToj. Therefore, Lu 2 Fe3Si5 



stands as the rare stoichiometric multigap superconduc- 
tor with the d-electrons in between the p-electron sys- 
tem (e.g., MgB 2 [1]) and the /-electron system such as 
PrOs4Sbi 2 3], providing an unique opportunity to study 
the MGSC in the moderately correlated electron system. 

The multigap superconductivity of Lu 2 Fe3Si5 is first 
observed by the specific heat measurement down to 0.3 
K under zero field [||. It is of interest to further elucidate 
the MGSC by performing the thermal conductivity mea- 
surements down to lower temperature in the vortex state, 
because the thermal conductivity is sensitive to the light 
carrier band which is expected to strongly affect on the 
low-energy QP excitations. Absence of the nuclear schot- 
tky contribution is another advantage of this technique 
in exploring the MGSC. 

In this paper, we report on the thermal transport mea- 
surements of single crystalline Lu 2 Fe3Sis down to T c /120. 
Our detailed results of the thermal conductivity in the 
mixed state confirm the MGSC in Lu 2 Fe3Sis. Moreover, 
from a comparative study with MgB 2 , Lu 2 Fe3Si5 reveals 
the significantly enhanced heat conduction as an indi- 
cation of the electronic interactions derived from the Fe 
3<i-electrons. 

Single crystals were grown by the floating-zone 



method [11] . One- heater-two-thermometer steady-state 
method was used to measure the thermal conductivity. 
The heat current q was aligned along the [001] direction, 
and the magnetic field is applied parallel to the ab plane. 
The thermal contacts with resistance of ~10 mil were 
made to the sample by using a spot welding technique. 

Figure [1] shows the temperature dependence of the 
thermal conductivity divided by temperature n(T) jT un- 
der zero field. An arrow denotes the superconducting 
transition temperature T c ~ 6 K. As clearly seen, n(T) JT 
shows a kink at T c followed by steep decrease with de- 
creasing temperature. On further cooling, a hump struc- 
ture appears around 3 K. The decrease of k(T)/T in the 
superconducting state attributes to the reduction of the 
QP density by opening the superconducting gaps. On 
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FIG. 1: (color online). Temperature dependence of the ther- 
mal conductivity divided by temperature k/T under zero 
field. The arrow denotes the superconducting temperature 
T c . Inset: k/T vs. T 2 plot under zero field. The solid line 
represents a fit to the data by k/T = kq/T + bT"' 1 . 



the other hand, the hump structure originates from an 
enhanced phonon mean free path due to the condensa- 
tion of electronic scattering centers in the superconduct- 
ing state. In fact, similar enhancement of the phononic 
conductance (called phonon peak) in the superconduct- 
ing state is also observed in various materials such as 
Nb [HI and Pb 0. The inset of Fig.[[]shows k/T vs. T 2 
plot. Since the measured thermal conductivity contains 
both electron and phonon contributions, it is necessary to 
separate to each part as k = K c + K p h- Below 0.15 K, k/T 
is well described by a relation of k/T = kq/T + bT 01-1 



with kq/T = 0.000 ± 0.005W/K 2 m, b = 0.11 W/K 2 8 m, 
and a ~ 1.8. Here, we should note that the conven- 
tional T 3 dependence of the phonon thermal conductiv- 
ity K P h dominated by the ballistic phonon scattering fails 
to reproduce our result. The lower power of K p h with 
a r~j 1.8 can be attributed to either the specular re- 
flection of the phonon or the phonon scattering off the 



electrons 14|, LL5J . For the electronic contribution, kq /T 



provides the electronic thermal conductivity in the zero 
temperature limit. The negligibly small kq/T clearly in- 
dicates that Lu2Fe3Si5 is a fully gapped superconductor 
consistent with the specific heat measurement 0. 

Figure [2] shows the field dependence of the thermal 
conductivity divided by temperature k(H)/T at several 
temperatures. For T > 1.0 K, n{H)/T takes a minimum 
at low fields and then increases with field. In general, 
the existence of the minimum of k(H)/T is explained as 
the result of the decrease of the phonon conductivity due 
to the vortex scattering concomitant with the increase 
of the electronic conductivity resultin g fro m the increase 
of the delocalized QP density 

0B1I21- Remarkably, 
for T < 1.0 K, k(H)/T shows rapid increase and reaches 
almost half of the normal state value K n /T already at 
low fields (poH < 1 T). Here, K n /T is estimated from 
the Wiedemann- Franz law via K n /T = Lq/po = 0.325 



W/K 2 m, where L = 2.44 x lO^Wft/K 2 is the Lorenz 
number, and pa = 7.5 pQcm is the normal-state resid- 
ual resistivity. It is of interest to compare our results 
with several superconductors. The inset of Fig. [2] depicts 
the normalized thermal conductivity k/k u at 0.1 K plot- 
ted against H/H c2 together with the data for Nb [l7j |. 
MgB 2 [HI, and UPt 3 [Hj]. Here, the phonon thermal 
conductivity k(H = 0) is subtracted from k/k u , and the 
upper critical field of poH C 2 — 6.4 T is obtained from 
Ref. [l9j . The field variation of k/k„ for Lu 2 Fe3Si5 is in 
dramatic contrast with the behavior of the fully gapped 
s-wave superconductor Nb [l7| , in which small magnetic 
fields hardly affect the low temperature thermal conduc- 
tivity. By contrast, in nodal superconductor UPt3, the 
delocalized QPs induced by the Doppler shift produces 
the remarkable field dependence of k/k u . The strongly 
enhanced thermal conductivity of Lu2Fe3Si5 is a clear in- 
dication of either a gap with nodes or nodeless multiple 
gaps with significantly different magnitude on distinct 
Fermi surfaces. However, the nodal gap is ruled out by 
the absence of the T-linear term of k(T). 

On the other hand, one immediately notices that 
Lu2FeaSi5 shares striking resemblance with MgB2, 
namely, a rapid increase of k/k u at low fields and a 
saturation behavior at high fields, although k/k u for 
Lu2Fe3Si5 shows even more pronounced increase and 
takes higher values. The field evolution of the thermal 
conductivity of MgB2 is well understood in terms of the 
multigap superconductivity with a small gap A s on one 
Fermi surface and a large gap A; on the other (index I 
and s represent large and small gaps, respectively) [lij ]. 
A consequence of the small gap is an existence of a "vir- 
tual upper critical field" H^ 2 , above which the overlap 
of huge vortex core provides a dramatic increase of the 
delocalized QPs contributed to the heat conduction. 

The characteristic field scale of H s c2 can be estimated 
from H s c2 /H c2 ~ (A s /A ; ) 2 (i>f, ; A<f>) 2 . On the other 
hand, the "normal state" contribution of the small gap 
band k s / n n is obtained from k s /k u = k s /(ki + k s ) = 
N s Vp s r s /(NiVp l Ti + N s v F s t s ) via a relation of k, cx 
NiV Fi Ti, where Ni, vp^, and Tj represent the normal 
state electronic densities of state, the Fermi velocity, and 
the scattering rate of each gap band (i = I, a). Here, 
suppose t s /t; is close to unity, k s /k u is simplified to 
~ 1/{1 + (Ni/N s )(vf,i/vf, s ) 2 }- Consequently, the two 
gap superconductivity is characterized by the three ra- 
tios A//A s , Ni/N s , and vf.i/vf.s- With the knowledge 
of A ; /A s ~ 4 Q and N^N,, v F ,i/v F , s ~ 1 0, we find 
H^ 2 /H C 2 ~ 1 and k s / K n ~ 0.5 for MgB2 in agreement 
with various estimations 

nana. 

Now, let us discuss the 
case of Lu 2 Fe 3 Si 5 . Since A//A s ~ 4 and N t /N s ~ 1 [| 
as same as MgB2, the k/k„ curve of L^FesSis is ex- 
pected to be scaled to that of MgB2 by tuning the ratio of 
vf,i/vf,s- We achieve excellent scaling with Vf,i/vf, 8 ~ 
0.8 as shown by th open circles in the inset of Fig. [21 
yielding H^ 2 /H C 2 ~ 0.04 and k s / 'n n ~ 0.6. The devia- 
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FIG. 2: (color online). Field dependence of the thermal con- 
ductivity divided by temperature k/T at several tempera- 
tures. Inset: The thermal conductivity normalized to its nor- 
mal state value /t//t n vs. H/H C 2 plo t at 0.1 K. For compari- 
son, the data for Nb [TJ, MgB 2 []]|, and UPt 3 are also 
shown. The open circles show a result of scaling described in 
the text. 



tion at high fields is due to the difference of H C 2 in these 
two systems. The inequality of Vf,i/vf,s indicates that 
the carrier mass m, which is inversely proportional to Vf 
(to oc 1/vp), is different in each band for Lu2Fe3Si5 in 
contrast to MgB2 [16J. This contradiction might origi- 
nate from the character of the dominant electrons con- 
tribute to the density of state at the Fermi surface; Fe 
3<i-electrons in Lu^FesSis and B 2p-electrons in MgB2, 
respectively. In addition, unequal weight of the d char- 
acter among the distinct Fermi surfaces yields the heavy 
and light carrier bands in Lu2Fe3Si5. Notably, a signa- 
ture of the electronic correlations possibly derived from 
the d-electrons can be found in the specific heat coeffi- 
cient j n — 23.7 mJ/molK 2 ||, which is about 10 times 
larger than that of MgB 2 (7„ = 2.62 mJ/molK 2 ) [Hj]. 
The rather small 7b a nd = 8.69 mJ/molK 2 obtained from 
the band calculations Q also indicates the presence of 
the electronic interactions. On the other hand, the small 
gap on the light carrier band provides dramatic enhance- 
ment of k(H)/T in the mixed state than that of MgB2 
because n is sensitive to the light carrier mass. 

In order to further clarify the multigap superconduc- 
tivity of Lu2Fe3Si5, we present the temperature depen- 
dence of k{T)/T in the vortex state (Fig. [3]). By apply- 
ing magnetic fields, we observe a pronounced increase 
of tt(T)/T at low temperatures corresponding to the 
fast growth of k(H)/T at low fields. Furthermore, an 
anomaly associated with a slight change of slope is found 
around 0.8 K at 0.1 T, and it shifts to - 0.3 K at 0.25 T 
as denoted by the arrows in Fig. [3] With further increas- 
ing fields, the anomaly becomes more pronounced, being 
a maximum around 0.2 K above 1.0 T. The emergence of 
the maximum can be attributed to either the phononic 
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FIG. 3: (color online). Temperature dependence of the ther- 
mal conductivity divided by temperature k/T under sev- 
eral magnetic fields. The solid and open arrows denote the 
anomalies associated with the small gap and T c at 5 T, re- 
spectively. The dashed line represents a fit to the data by 
k/T = k /T + bT a ~ 1 . Inset: The calculated k/T curves for 
the corresponding fields within the two-gap model described 
in the text. 



or the electronic contribution. However, we immediately 
exclude the phononic origin because the appearance of 
phonon peaks twice below T c is highly unlikely. More- 
over, an increase of the phonon mean-free path up to 
10 folds of the zero-field value, for example at 5 T, is 
improbable even if the low-temperature thermal conduc- 
tivity is dominated by the specular reflection. Thus, we 
turn to the electronic origin raising the following possi- 
bilities; (1) an enhancement of the QP mean free path l e , 
(2) a magnetic anomaly associated with a magnetic order, 
and (3) an increase of delocalized QPs excited above the 
field-suppressed small gap A s . For the possibility of (1), 
the enhancement of l e upon entering the superconducting 
state usually occurs as a result of the strong inelastic scat- 
tering. This is because normal state l e suppressed by the 
inelastic scattering starts to increase in the superconduct- 
ing state due to the condensation of electronic scattering 
centers. In practice, this behavior has been observed in 
CeCoIns [22| and YBCO (23|, in which a source of the 
inelastic scattering is attributed to the magnetic fluctua- 
tions. However, no signature of such magnetic fluctuation 
has been indicated in Lu2Fe3Sis [24(. The magnetic ori- 
gin of the maximum is also ruled out because Fe atoms 
carry no magnetic moment [24[ . Therefore, in the follow- 
ing we discuss the remaining possibility (3) by analyzing 
k(T)/T within the framework of the two-gap model. 

In the superconducting state, the electronic thermal 
conductivity tt e is expressed as, 



_ 2F 1 (-y) + 2yln(l + e~») + y 2 (l + e^)" 1 

Ke/Kn- 2Fi(Q) , (1) 
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where y = A(t)/kB,T, A(t) being the half of the energy 
gap and t = T/T c [25[. The temperature dependence 
of A(i) is given by standard gap interpolation formula 
A(t) = A tanh(2.2 A /T7* - 1)> where A 

o is an energy 

gap at absolute zero. The term of y) is given by the 



expression F n (—y) — J °° z n (l 



i z+v )~ 1 dz. Following 
the approach used for the specific heat Q , we generalize 
Eq.([T]) to the two-gap model as n e = n s n etS + (1 — n s )K et i, 
where K e ^ (i = I, s) is the large and small gap bands ther- 
mal conductivity, respectively, and n s is a weight for the 
small gap band. In the following analysis, n s is fixed to be 
0.6 regardless with the applied field, which is determined 
from the scaling. To obtain c„, we estimate the phonon 
thermal conductivity from the minima of k(H)/T shown 
in Fig. [2] assuming that the minima represent at most the 
maximum value of the phonon contribution [l6| . 

The solid lines shown in the inset of Fig. [3] repre- 
sent the calculated results of the total thermal con- 
ductivity k/T = K e /T + Kph/T based on the two-gap 
model at each field. Here, we assume the temperature 
dependence of the phonon thermal conductivity to be 
Kph/T = 0.11 x T ' 8 as determined by the zero-field k/T 
(the dashed line in Fig. [3]). Moreover, at zero field, we use 
the gap values of 2A 0t i/k B T c = 4.4 and 2A , s //c B T'c = 1-1 
obtained from the specific heat measurement |8(]. The 
field variation of Ai(H) is assumed to follow the mean- 
field description Ai(H) = Ao^-^/l — H/H C 2 (the solid 
line in Fig.0]). It should be emphasized that the two-gap 
model well reproduces the experimental results especially 
for the maximum of k(T)/T only by tuning the small gap 
A s /fcBT c . The result further supports the existence of 
the small gap in Lu2Fe3Si5. Furthermore, the maximum 
most likely originates from an increase of the delocalized 
QPs excited over the small gap. On the other hand, there 
exists a discrepancy between the experiment and the cal- 
culation below the maximum temperature; the experi- 
mental results show gradual decrease while the calculated 
k/T rapidly drops to zero. One possible interpretation 
is that the system behaves like a dirty superconductor 
(£ > l e ) due to the presence of n H^ 2 that gives rise 
to a large coherence length £. It has been argued that 
the superconductors in the dirty regime shows a rapid 
growth of QP density in the mixed state in comparison 
with those in the clean limit (£ <^ e ) 

The parameter of A s /k B obtained from the analy- 
sis and A;/fce are plotted against the magnetic field in 
Fig. 01 It is clearly seen that A s /fce sharply decreases at 
low field fJ,oH s < 0.25 T, while Ai/k-Q shows monotonous 
decrease. Note that the characteristic field scale of hqH s 
is close to the virtual upper critical field HoH^ 2 ~ 0.26 
T. Interestingly, a similar field-induced suppression of A s 
below HoHc2 ls a ^ so observed in MgB2 as demonstrated in 
the inset of Fig.Hl in which t he g ap values are determined 
by the point-contact study [2j|. For the case of MgB 2 , 
this behavior is understood as a consequence of a weak 
interband pairing interaction 27|, |28[ . From the analogy 
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FIG. 4: (color online). Field dependence of the supercon- 
ducting gaps A/ks- A(//cb and As/ks represent the large 
and small gaps, respectively. The solid line denotes the cal- 
culated result of Aj/fce based on the mean-field description. 
The dashed line is a guide to the eye. Inset: The normal- 
ized small gap A a /A S (H = 0) vs. H/H^ 2 for L^FesSis and 
MgB 2 [13]. 



of MgB2, the existence of the multiple bands having the 
different dimensionality and the weak interband interac- 
tion could be a source of the MGSC in Lu2Fe 3 Sis. 

In summary, our thermal conductivity measurements 
clarify the multigap superconductivity with a fully 
gapped excitation spectrum in L^FeaSis. The weak cou- 
pling between the distinct gap bands is thought to be 
the origin of the multigap superconductivity. In contrast 
to MgB2, the dramatic enhancement of the quasiparti- 
cle heat conduction in the mixed state implies the pres- 
ence of the electronic correlations derived from the Fe 
3d-electrons. Our findings shed new light on the multi- 
gap superconductivity in the correlated electron systems. 
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